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Abstract 

Quark models with four-fermion interaction including derivatives of fields in the 
strong coupling regime are used to implement composite-Higgs extensions of the 
Standard Model. In this approach the dynamical breaking of chiral symmetry occurs 
in two (or more) channels (near polycritical values for coupling constants), giving 
rise to two (or more) composite Higgs doublets. Two types of models are built 
for which Flavour Changing Neutral Currents (FCNC) are naturally suppressed. 
In the first Model I the second Higgs doublet is regarded as a radial excitation 
of the first one. In the second Model II the quasilocal Yukawa interaction with 
Higgs doublets reduces at low energies to a conventional local one where each Higgs 
doublet couples to a definite charge current and its v.e.v. brings the mass either 
to up- or to down-components of fermion doublets. For the special configuration 
of four-fermion coupling constants the dynamical CP-violation in the Higgs sector 
appears as a result of complexity of v.e..v. for Higgs doublets. 
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1. Introduction 



The fundamental particles of the Standard Model (SM) of electroweak interactions, 
leptons, quarks and gauge bosons, acquire masses through the interaction with a scalar 
field (Higgs boson). The mass generation is mediated by the Higgs mechanism which rests 
on the Electroweak Symmetry Breaking (EWSB). To accommodate the well- established 
electromagnetic and weak phenomena, the Higgs mechanism requires the existence of at 
least one iso-doublet, complex scalar field. After absorbing three Goldstone modes to 
build massive states of W^,Z bosons, one degree of freedom remains, corresponding to 
a real scalar particle (Higgs boson). If SM holds valid as a weak-coupling theory till 
very high energies then this particle cannot have a mass heavier than few hundreds of 
GeV. Thus the search for the Higgs boson is one of the fundamental quests for testing 
the minimal SM. Current estimations based on the different theoretical requirements and 
experimental implications JI], ^ give the SM Higgs mass in the "intermediate mass" window 
65 < Mh < 200 GeV for a top quark mass value of about 175 GeV [^. Despite of the 
recent successes of the SM in its excellent agreement with the precision measurements 
at present energies [^, it is generally believed that the SM is not the final theory of 
elementary particle interactions. 

There are many extensions of the SM which lead to the enlargement of the Higgs sector 
of the SM. For instance, the Minimal Supersymmetric Standard Model (MSSM) entails 
two elementary Higgs doublets at low energies, the Two-Higgs-Doublet Model (2HDM) 
contains two complex SU{2)l - doublet scalar fields with hypercharge Y = ±1 to couple 
the up-type / down-type right-handed quarks to its Higgs doublet. The search for relations 
between the many Higgs-field dynamics and the masses of t-quark and Higgs boson give 
the selection rule for a particular model beyond the SM as well as for its acceptable 
parameters 0, 0, f^, In more complicated theories ( see and references therein) 
such as SUSY SM ones, Eq ones, or Left-Right symmetric ones [0, several neutral 
scalars, charged scalars and even double-charged scalars are required in order to give all 
amplitudes acceptable high-energy behavior [jl], 0. 

However there exists an alternative possibility [|ri[] to restrict the number of elementary 
particles to the observable fermion and vector-boson sector with generation of scalar Higgs 
particles due to attractive self-fermion interaction. Namely, the quark self-interaction may 
be responsible for the production of quark- ant iquark bound states which are identified 
as composite Higgs particles. The idea that the Higgs boson could be a bound state of 
heavy quark pairs has been developed and worked out in a series of papers being 
motivated by the earlier work of Nambu and Jona- Lasinio (N JL) . In particular, for 
t- quarks, it is provided by the Top-Mode Standard Model (TSM) Lagrangian, known also 
as the Bardeen-Hill-Lindner (BHL) Lagrangian |11| , , . The possibility that multi- 
ple four-fermion interactions (for three and a heavy fourth generations) are important in 
EWDSB, leading to an effective 2HDM at low energies, has been investigated in In 
this model Higgs boson induced Flavour Changing Neutral Currents (FCNS's) are natu- 
rally suppressed |T^. Some recent theoretical aspects and questions of it- condensation 



frameworks one can find in the review of |17|. In these scenarios the heavy top mass is 



explained by the "top-condensation" where new strong forces lead to the formation of tt 
bound states and the EWSB. In a minimal version of quark models the top-condensation 
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was triggered by a local four-fermion interaction. 

The main goal of this paper is to give the description of the design of the Quasilocal 
Quark Models of type I and type II which provide two composite Higgs doublets and 
in principle satisfy phenomenological restrictions on FCNC suppression. The particu- 
lar, sample choice of formfactors and bare Yukawa coupling constants is made to obtain 
estimations for typical mass spectra in Model I and Model II. 

We propose the quark models with Quasilocal four-fermion interaction |TB[ where the 



derivatives of fermion fields are included into vertices to influence on the formation of the 
second Higgs doublet. Such extensions of the Higgs sector lead to a broad spectrum of 
excited bound states, moreover they may be viewed as more natural than other, above 
mentioned extensions since the particles involved in EWSB form only a ground state spec- 
trum generic for SM. In these quasilocal NJL-like quark models (QNJLM) the symmetries 
do not forbid further higher dimensional vertices and one should expect that the ground 
states could be accompanied by (radial) excitations with identical quantum numbers but 



much higher masses |T9|, ^ ^ p^] . 

Thus, from the viewpoint of the 2HDM SM, the QNJLM are attractive because: i) 
it is an extension of the minimal TSM which adds new phenomena (e.g. a broad mass 
spectrum of bound states including charged Higgs bosons); ii) it is a minimal extension 
in that it adds the fewest new arbitrary constants; iii) it easily satisfies theoretical con- 
straints on p ~ 1 and the tree-level FCNC's suppression in accordance with the 
experimental evidence; iv) such a Higgs structure is required in order to build a model 
with the CP- violation ||23| because the one-Higgs doublet interaction does not provide 
any effect of dynamical CP- violation. We shall show in a toy model with quasilocal 
four-fermion interaction how P-parity breaks down dynamically for the special choice of 
coupling constants 



This article is organized as follows: Section 2 contains the simplest Gross-Neveu model 
which reminds how the Dynamical Chiral Symmetry Breaking ( DCSB) arises in the 
scalar channel due to strong interaction. In Sec. 3, we formulate the main rules for 
the construction of the QNJLM which admits the polycritical regime. Here the effective 
potential and the mass spectrum for composite scalar and pseudoscalar states are derived 
for them. For more evidence, in Subsec. 3.3, we investigate the two- channel QNJLM in 
the large-log approximation. In the vicinity of the tricritical point all possible solutions 
are analyzed. It turns out that there exist three phases with different correlation lengths 
in the scalar channel. Moreover the special phase of dynamical P-parity breaking is 
found. In Sec. 4 and 5 two types of models are built for which FCNC suppression may be 
naturally implemented. In Sec. 4, the first extension of the SM composite two-Higgs bosons 
for QNJLM (2HQM) is proposed where the second Higgs doublet is regarded as a radial 
excitation of the first one. The second model is constructed in Sec. 5 so that the quasilocal 
Yukawa interaction with Higgs doublets reduces at low energies to a conventional local 
Yukawa vertex where each Higgs doublet couples to a definite charge current and its 
v.e.v. brings the mass either to up-or to down-components of fermion doublets. In this 
version the top and bottom masses are explained by "top-, bottom- condensations". On 
the base of the effective potential for the Model II and the positivity of its variation the 
mass spectrum for composite states is investigated. It is interesting that for the special 
configuration of coupling constants the dynamical CP- violation appears in the Higgs 
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sector. In the Summary we discuss the obtained resuhs and a possibihty to use them in 
different aspects of high energy physics. The Appendix contains the calculation of the 
matrices of the second variations for composite two-Higgs bosons in Model I ,11 and the 
effective potential of Model II for the special choice of quasilocal formfactors. 



2. DCSB in Models with Four-fermion Interaction 
and the Critical Regime 

Let us remind how the Dynamical Chiral Symmetry Breaking (DCSB) arises in a 
model with local 4-fermion interaction due to strong attraction in the scalar channel. The 
simplest, Gross-Neveu (GN) model retaining the scalar channel only can be presented by 
the Lagrangian density, in two forms (the Euclidean-space formulation is taken here), 

C = ql/)q + -r^imf = q{V> + ^0(a;))g +^02(x), (1) 

where p = i7^<9^ and q = (qi) stands for color fermion fields with Nc components. For 
the time being we take the number of flavours Np = 1 and the current quark mass rriq = 0. 
In Eq.(|l]) the scalar auxiliary field 0(x) (a prototype of the Higgs field) is introduced in 
order to describe the dynamical symmetry breaking phenomenon in the large- A^^c limit. 

This model is implemented by an 0(4)-symmetric momentum cutoff A for the fermion 
energy spectrum. For a quark model the cutoff A can be thought of as a separation 
scale which appears when evaluating the SM low-energy effective action from a more 
fundamental theory. The regularized effective action Seff for auxiliary field, 

Z^(0) = expi-S^jf) = ^exp(- |d4x/:(0(x)))^_^, (2) 

possesses the mean-field extremum on constant configurations =< >= rud = const. 
The relevant effective potential V^jf can be obtained by integration over fermions, 

(3) 

where the constant /i is a normalization scale for quark fields. Its extrema can be derived 
from the mass-gap equation. 

The main contribution into Eq.(^ is given by a tadpole term in the fermion loop which 
is related to a vacuum expectation value (v.e.v.) of the scalar fermion density, 

, , Svr^A^ 47r\ , 
/?(0)= 0— ^ + ^— (gg). (5) 
9 'c 
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The cutoff independence is realized with aid of fine-tuning, 871^/(7^ ~ 1 — 0(1/A^). In 
the language of the theory of critical phenomena it is equivalent to developing our model 
around a critical or scaling point where the quantum system undergoes the second-order 
phase transition. By definition the critical coupling constant is g1^^^ = Svr^. When g"^ < 
g1^^^ the only solution of mass-gap Eq.(^ is = 0, while for g"^ > g1^^^ there exists another 
nontrivial solution for dynamical mass 7^ which brings the true minimum for Ve//. 
Meanwhile the symmetric solution = does not provide then a minimum anymore but 
realizes a maximum. 

The fine-tuning states that the strong A^-dependence should be compensated by the 
corresponding term in the coupling constant, 

Its practical meaning is evident, namely, one produces a mass scale for physical states 
which is much less than the cutoff scale governing large radiative corrections. The devi- 
ation scale m-Q << A^ determines the physical mass of scalar meson. Indeed its kinetic 
term can be obtained from the second variation of S^jf by calculating the 1-fermion loop 
diagram (see App. Fig.l), 

1 /• rfV 7 



Seff ^ 5e//(0 = ^d) + </)(-p)r(p) (f){p), 

(f) = nid + (f>, (7) 
where the inverse propagator of scalar field reads: 

g^ Jk<A {27ip L 

= K + pW) + o(^ 

in the chirally invariant regularization of the fermion loop. The scalar meson mass is given 
by the remarkable Nambu relation m<^ ~ 2md and the formfactor I{p) is determined by 
the relation, 

OAT f ^'^^ ^ ^ 

I[p) = '^N.j^^^ ^2^^^ ^ ^^^2 + m'i{k- \pY + ml' 

In order that the physical mass parameters were insensitive to A, i.e. d\md = 0, the 
scale mo should be weakly dependent, ml ~ ln(A2/m^), on the cutoff A. 
What have we learned from the above model? 

(i) The cut-off theory can be used for processes involving momenta p much less than A 
for the purposes of discarding high-energy states from the theory. 

(ii) The mass scale of meson states is assumed to be much less than A which is imple- 
mented in the vicinity of critical values of coupling constants. 

(iii) As a result of DCSB in these models only one type of scalar mesons (i.e. eventually 
one Higgs doublet) is created in the large- iVc approach. 
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(iv) In such models the (radial) excitations of composite meson states are not present in 
the la.Tge-Nc approach. 

Meantime the conventional quark models with local four-fermion interaction may not 
represent a consistent part of the Beyond-Standard Model (BSM) effective action and 
conceivably they shall be extended with inclusion of higher dimensional vertices which 
are not forbidden by symmetries and induce the appearance of a reach spectrum of excited 
composite meson states. 



3. Quasilocal Quark Models and Polycritical Regime 

3.1. Dominant higher-dimensional vertices in DCSB phase 

In order to involve in the theory the effects of the discarded states at scales of order A 
it is needed to adjust the existing couplings constants in the Lagrangian and to add new, 
quasilocal, non-renormalizable interactions (vertices). These vertices are polynomial in 
the fields and derivatives of the fields and only a finite number of interactions is required 
when working to a particular order in ^, where p is a typical momentum in whatever 
process is under study. 

We examine the DCSB patterns in the mean-field approach (large-A^^c limit) and es- 
timate the vertices with any number of fermion legs and derivatives. The main rule to 
select out relevant vertices is derived from the requirement of insensitivity in respect to 



the separation scale A following the conception of low-energy effective action [|T8 
We assume that: 

(i) A^-order contributions from different vertices are dominant in creating the DCSB- 
critical surface that is provided by cancellation of all contributions of A^-order and 
defines the polycritical regime; 

(ii) A°-order contributions from vertices assemble in the mean-field action to supply 
fermions with dynamical mass rrid << A which establishes the low-energy physical 
scale; 

(iii) respectively A^^ (etc.)-order contributions are irrelevant at energies much lower than 
A and so may be dropped from the theory if such accuracy is unnecessary. 

In the large-A^^c approach the following approximation for v.e.v. of fermion operators is 
valid, 

= + 0(l/iV,)), (10) 

where any number of derivatives can be inserted between antifermion and fermion oper- 
ators. 

V.e.v. of a bilinear operator is estimated in the assumption that quarks obtain a 
dynamical mass. Namely, 

^ I ^ tr-f—. N^nidA^. (11) 



A2" J (27r)4 fS + imd 

|p|<A 
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One can see that the vertices with derivatives in many-fermion interaction are not sup- 
pressed and play equal role in the mass-gap equation. 

We omit the full classification of effective vertices relevant in the mass-gap Eq. (see 
IS]) and report only the minimal structure of the QNJLM which admits the polycritical 
regime, 

1 ' _ ( -d'^\ - /-<9^\ 

C = qpq+ ^rnn qufn j ' ^^fm j^^?' (^2) 

where a^n is a hermitian matrix of coupling constants without zero eigenvalues and it is 
taken to be real symmetric one in order that the interaction did not break the CP-parity 
explicitely. Chiral fermion fields are given by qL(R) = 1/2(1 ± Ts)^'- We define the vertex 
formfactors to be polynomials of derivatives, 

Ur)=Ef^^r\ (13) 

i=0 

to have quasilocal interactions. The variable r is related to derivatives, r —d^/h?. 
We adopt the following rule for derivative action which provide the hermiticity of fermion 
currents: 

Qt^Q = 7? — I — (14) 



^A2^"4"V A 

Besides let us regularize the interaction vertices with the help of a momentum cutoff, 

qq > qe{K^ + d'^)q. (15) 

Without loss of generality one can choose formfactors /^(t) being orthogonal polynomials 
on the unit interval, 

1 

dTfm{j)fn{T) = 5mn- (16) 

Let us now introduce the appropriate set of auxiliary fields 4>n{x) ~ const and develop 
the mean-field approach, 

£(0) = q(p + zM{<P)PL + zM^i4>)PR)q + N,A' E 0™a™„0n. (17) 

^ ^ m,n=l 

The dynamical mass functional is a linear combination of formfactors, 

M(0) = E0.(a;)/„f^). (18) 



n=l 



In accordance with Eq.(p!^ the differential operator M{(f)) is understood as a Weyl ordered 
or fully antisymmetrized product of functions (pn and derivatives. Thereby we come to a 
model with I channels. When integrating out the fermion fields one obtains the effective 
action of 0*, - fields. The effective potential Veff is proved to be a functional depending 
on the dynamical mass functional M{(j)*, 0) and proportional to Nc that allows us to use 
the saddle point approximation for Nc » 1. 



7 



3.2. Effective potential and equations on the mass spectrum for 
QNJL Model 

The effective potential for auxiliary fields can be derived with the momentum cutoff 
regularization by averaging over quark fields: 



' N } I 

iV,A2 ^ 0;,a->„ - dTT\nil + 

m,n=l Q V 



|M(r)p 
AV 



/V r ' 1 / 1 



2' - V |MoP 2 



VA2 



--J^{\M{r)\'-\MX) + 



(19) 



herein Mq = M(0). The last approximation is valid in such a strong coupling regime 
where the dynamical mass Mq << A. This regime is of our main interest and it is 
realized in the vicinity of a (poly) critical surface. The critical values of coupling constants, 
^mn — ^mn/Svr^, are found from the cancellation of quadratic divergences. In this paper 
we study the critical regime in all / channels. The vicinity of this polycritical point is 
described by the following parametrization: 



A 

A2 

The generalized mass gap equations, 



7i^a-i ~ S^n + -i^, I A,,| « A'. (20) 



deliver the extremum to the effective potential which may cause the DCSB if it is an 
absolute minimum. They read: 



^ r (if 
Y.Arnn<Pn = J ^j^{\M{t)\' M{r) 



dr 

n=l ~^ 



A 



2 



/™(0)|Mo|^Moln— 2 



+ ' 



j -^{\M{T)\'M{T)Ur) - rn'MoUm. (22) 



It can be seen from the first relation that, 

^ CA^„0„ = I dr > 0, (23) 

m,n=l g T -\- -j^ 

which means that for the existence of a non-trivial dynamical mass it is necessary to have 
at least one positive eigenvalue of the matrix A^n- However not all the solutions provide 
a minimum (see, the analysis of two-channel models in p4|, ||19|| ). 
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The true minimum is derived from the positivity of the second variation of the effective 
action around a solution of the mass-gap equation, 



<Pm=<<Pm> +(ym{x) + i-nm{x). 



(24) 



This variation reads 



+2 (vr, (A"V + 5"")(t) + (tt, (A'^V + B^^)^ 



(25) 



where two symmetric matrices - for the kinetic term A = (v4^„) , i,j = (o", tt) and for 
the constant, momentum independent part, B = (-B,^„) - have been introduced. 

The positivity of the second variation corresponds to the formation of physical mass 
spectrum for composite scalar and pseudoscalar states which can be found from zeroes of 
the second variation determinant at the Minkovski momenta (p^ < 0), 



det(i/ + 5) = 0. 
Matrix elements of B are given by the following relations: 



(26) 



B" 



Bl 



6 



dr 

T 



ReMff^{r)Ur) 



+ MoV™(0)/„(0)(61n^-4 



M2/m(0)/„(0) 
- 2A 



1 



+ 2 



dr 



2 — 



ReM)2/™(r)/„(r) - M2/™(0)/„(0) 



+ 2MoV,n(0)/„(0)ln— 2 -2A, 

Mq 



(27) 



(2^ 



/d-T 
— (ReM)( ImM)/^(r)A(r), 



(29) 



where the terms of 1/ A^-order are neglected. 

When exploiting the mass-gap equation (^) one can prove that the matrix 13 has 
always a zero eigenvalue related to the eigenvector =< vr^ > —i- < cr^ >. It corresponds 
to the arising of the Goldstone mode (the massless Goldstone bosons). 
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The kinetic energy matrix A turns out to be block-diagonal 

1 



mn "5 mn mn 



1 + 



In^ 



(30) 



A. 



aa 

mn 



+ 



/^(0)/„(0) \\^ — ^0{\] 
[UT)fnir) - /„(0)/„(0)] 



dr 



+ 



A2 



(31) 



herein we have displayed the leading terms only in the large-log approximation. The more 
detailed expression can be found in |T^ . 



3.3. Quasilocal Two-Channel quark models and possibility of 
dynamical breaking of P-parity 

For the further investigation of composite Higgs extensions of the Standard Model let us 
consider the Quasilocal Two-Channel quark model in a tricritical point [|TP|]. We set 
m,n = 1, 2 in (12)-(20) and retain only the lowest derivatives in the potential, with /i = 1, 
/2 = a/3(1 — 2r). The dynamical mass function is thereby, M(0) = 0i + 02\/3(l — 2r). As 
(pj are complex functions, M{(j)) is complex too. However, with the global chiral rotation 
M[(j)) —>■ M(0)e*'^, uj = const it is always possible to implement Im < Mq >,/,= and 
we can choose the following parameterization: 



>i + ip, 



(pi = Re0i 



(32) 



The equations (|21| ) for the Two-Channel model read: 

A2 



Aii0i + Ai202 
C?10i - ^202 

piVSAu - A12) 



In ^ - 6 V3</.?02 - 18ct>i4>l - 8^30: 

4 



Mi 



2^301(02 



20: 



2), 



2pV3{4>l + <Pl + ^p2^ 



where 



di = VSAn - A 



12) 



dy. 



-V3A 



21 



A 



22- 



(33) 



(34) 



We analyze the equations (|33D near a polycritical point, |Ajj| ~ ^ A^, in the large- log 



approximation (In ^ ^ Inln ^). It gives rise to a set of solutions. 
For p = all the solutions are divided into the following classes: 



a) Gross-Neveu-like solutions 



lGN 



are: 



dl det A 



{V^d^ + d2r\n^ 



1 + 



In^ 



di 



(35) 
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These solutions deliver minima to the potential when ^/3dl + 6,2 < 0, with one eigenvalue 
of the matrix A being in the over-critical regime and the other one in the sub-critical, 
b) Abnormal solutions are: 



1 + 



^ -— (36) 



they correspond to the suppression of the large log-terms in Eqs.(33) of motion and give 
minima to the potential, when \pid\ -|- ci2 > 0, \/l>d\ — 2d2 7^ ( either both eigenvalues 
of A are positive, or one is positive and the other one is negative ). 

c) On the planes \/?)di + ^2 = and y^rfi — 2d2 = there appear special solutions 
with different, peculiar asymptotics |^, |T^. 



d) In general, in the models with more than one channel complex solutions are allowed, 
and the imaginary parts of all the variables (pj cannot be removed simultaneously by a 
global chiral rotation. However the complex solutions (p 7^ 0) minimize the effective 
action only (!) for the narrow domain in the vicinity of the plane \/3di — 2^2 = 0. Their 
asymptotic expressions are: 



l2 



di + 4Ai2 



-Vs^i, (37) 



16v/3(ln^-3) 

and the dynamical mass is = 40^. The axial part of the mass function looks as follows: 



2 di\/3 3 2 , !2\ di\^ 



1 + 



In^ 



(38) 



In each of the phase space domains mentioned above one finds four common boson 
states — two scalar and two pseudoscalar — for real 0^, and, in general, — for complex 
0j, three states with mixed P-parity and the pseudoscalar one with zero mass, the latest 
is in accordance to the Goldstone theorem. 

The mass spectrum of related bosonic states (collective excitations) is determined 
by zero-modes of the matrix of second variations of the effective potential (25) and respec- 
tively by Eqs. (26) - (31). Taking into account the conditions necessary for a minimum 
of the potential, we find the solutions at —m^ = < 0, giving physical values of particle 
masses. 

In the case of p = 0: 

a) NJL-like mass spectrum: 

V^c/i + d2 



2 2 2 V oui a2 
m„ = U m, ^ m, ^ 



< ^ 4m^, (39) 

in this domain the radial excitation states are heavier than the lightest scalar meson by 
a factor of logarithm. 
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b) For the Abnormal solutions we have: 



m^=0 ^^^iflV'M*!:^^— (40) 



2 

o 



When comparing (39) and (^Op we find the scalar channel correlation length to be different 
for each phase, that corresponds to the tricritical point conditions. 

c) For the special real solutions the relations between scalar and pseudoscalar meson 



masses are different from (p9D,(^0D (see [p^ ],|T9 



d) Mass Spectrum in the P-parity Breaking Phase ( p 7^ 0). One can see from 
that in the large-log approximation the axial dynamical mass (the imaginary part of 
M(0)) dominates. It leads to appearance of a massless boson in the scalar channel in 
accordance to the Goldstone theorem. Conventionally, the massless boson is related to be 
a pseudoscalar meson corresponding to the generation of a real dynamical mass. In order 
to fit it we make a global chiral rotation of fermionic fields q exp(i757r/4)g accompanied 
by corresponding rotation of the bosonic variables (f)j i(f)j: 

01 = -p, 02 = i(p2 + (41) 

The classification of states given by the P-parity quantum number is relevant only in the 
large-log approximation, when: 

^^^ = 0[ I , (42) 

next-to-leading logarithmic effects are of no importance and one can neglect the mixing 
of states with different P-parity. Then the mass spectrum of mesons is: 



di + 4Ai2 



mf = 0, ml^ "^^^""22^ ~ 1601 = 



V3d, m' ^ ^(^^ + ^''^ m 



The ratio of m2 and m^ does not depend on the logarithm, so both the masses are 
comparable. On the other hand, in the models with a finite momentum cut-off, when the 
effects of order of 1/ In ^ make sense, the dynamical P-parity breaking is induced, since 
7^ 0. This phenomenon of dynamical P-parity breaking can be used in extensions of 



the Standard Model |TT[ where several Higgs bosons are composite ones. 

Thus we conclude that the models with polycritical (tricritical) points are drastically 
different from the local NJL models in the variety of the physical phenomena in the DCSB. 
Explorations of such QNJLM in extensions of the SM are pretty well motivated as the 
underlying dynamics responsible for the top quark condensate should most likely lead 
to a broad spectrum of excited states, just like the hadron dynamics with QCD as an 
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underlying force. Moreover the QNJLM may even be viewed as more natural than the 
extensions to more generations, more elementary Higgses, or to SUSY in the SM since 
the particles involved in DCSB (with masses of order of the electroweak scale) belong in 
this context only to the ground state spectrum. In the following sections we shall present 
an extension of the SM with two-Higgs bosons of the QNJLM where one of the Higgses 
is a radial excitation of another one. 



4. Higgs Bosons as Radial Excitations - Model I 
4.1. Effective potential in Model I 

Let us construct now the two-fiavor quark models with quasilocal interaction in which the 
t- and 6-quarks are involved in the DCSB. In accordance with the SM, the left components 
of both quarks form a doublet: 

'tL 



U = [ J , (44) 

which transforms under SU{2)l group as a fundamental representation while the right 
components tji, hn are singlets. 

The Model I which is to satisfy the FCNC suppression has the following Lagrangian: 



an (gt,kJt,k + 9b,kJb,k) {gb,iJb,i - gt,iJt,i) ■ (45) 



Here we have introduced the denotations for doublets of fermion currents: 

Jt,k = tRft,k [~j2jlL, Jb,k = bRfb,k I 1 Ql, (46) 

and the tilde in Jt± and Jb^k marks charge conjugated quark currents, roteted with r2 
Pauli matrix 

Jt,k = 1^-2 Jt^k^ Jb,k = ir2Jlk (47) 

The subscripts t, b indicate right components of t and b quarks in the currents, the index 
k enumerates the formfactors: 



As the spinor indices are contracted to each other in (|46D, Jt^k transforms as a doublet 
under SU{2)l. T2 is a Pauli matrix in the adjoint representation of the group SU{2)l. 
Coupling constants of the four-fermion interaction are represented by 2 x 2 matrix aui and 
contributed also from the Yukawa constants gk,i,gb,k- 
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The Lagrangian density of the Model I (|45|) to describe the dynamics of composite 
Higgs bosons can be obtained by means of introduction of auxihary bosonic variables and 
by integrating out fermionic degrees of freedom. According to this scheme, we define two 
scalar 5'f7(2)i-isodoublets: 



$1 



and their charge conjugates: 

$1 



$5 



$2 



021 
>22, 



^22 
021, 



(49) 



(50) 



In terms of auxiliary fields, the Lagrangian ( ^5]) can be rewritten in the following way: 

iVA2 2 , 2 

k=l 



J 



L 



kin 



E + ^ E [9t,k^Jt,k + gb,k'^Ub,k 



k,l=l 



h.c. 



(51) 



The integrating out of fermionic degrees of freedom will produce the effective action for 
Higgs bosons of which we shall keep only the kinetic term and the effective potential 
consisting of two- and four-particles vertices. The omitted terms are supposedly small, 
being proportional to inverse powers of a large scale factor A. The Yukawa constants are 
chosen of the form 

gt,k = 1; 9b,k = g (52) 

for k = (1,2). The first constant is set to unit because the fields $i and $2 can always 
be multiplied by an arbitrary factor which is absorbed by four fermion coupling constants 
through redefinition of polycritical coupling constants. The remaining constant g induces 
thereby the quark mass ratio mh/mt. The effective potential for the Model I has the 
following form: 



eff 



^ f- E A.. + 8($t$0^ fin 

Svr^ V \ 4($l$i) 



+ 



1 



+ 



k,l=l 
4^^t^ \2 



8(?"($1<1'1 



In 



A^ 



9(l + ^^)/^t^ ^2 3(1+^4) 



4 



(<I>I<l>l)($^$2) + 



3(1 + 9^) /^t^ ^/^t^ N 3(1 + g^) _x ^2 3(1 + g^) ..^ 

- / ^ ($I$2)($^^i) + ' (^1^2)' + ^-—^m^iY 

4 00 



4 



($I$2) + {^1^1 
$I$2) + ($^$1) 



o 



InA^ 



(53) 



where the bilinear "mass" term is in general non-diagonal and represented by the real, 
symmetric 2x2 matrix Am- 
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We assume the electric charge stabihty of vacuum or, in other words, that only neutral 
components of both Higgs doublets may have nonzero v.e.v. Hence, one can deal with only 
neutral components of the Higgs doublets in the effective action for studying DCSB. This 
part of the Higgs sector can be investigated separately as a model where two singlets (not 
doublets) appear as composite Higgs bosons. For this purpose, we use the Quasilocal Two- 
Channel model which we have already developed in Sec. 3 for the case of one-flavour [|19|. 

Following the definitions made in [^|, we relate the fields 0i, 02 and p to the neutral 
components of Higgs doublets: 



011 = 0i; 022 = 02 + ip- (54) 

The condition of minimum of the potential (^) with the charged components of Higgs 
doublets put to zero values: 0i2 = 02i = 0, brings the mass-gap equations for them: 



Aii0i + Ai202 



160? In ^ + 16/0? In - {I + g^' 



159 



15^3 



'102 



+ ^0102 + ^02 + ^(301 + V^02 



Ai20i + A 



22V^2 



:i+^?^)[-^0? + |0?02 + ^0i0^ + 



9.3 



^(v^0i + 902) 



4pA22 = p(l + g^) (301 + 2^30102 + 90^ + 9p^ 



Let us consider the equations (|55D-(p7D for two cases: 1) p = and 2) p 7^ 0. 

When p = 0, assuming that 0i 7^ and 02 7^ 0, we rewrite the equations (|55|) -l 
the following way: 



(55) 

(56) 
(57) 

in 



All 



?ln 



^^,+160?,^ln^^,^, 

\/3 0i 



+ 402 



A22 = (1 + 

9.2 



+ {i + g' 



159 



15^3 



>W2 



+ 



4 01 
^\ [9 

5v^0? 



A12— , 



(5^ 



1 + 



3^3 



+ 



4 



A 

^12 — 



The solution of the mass-gap equation of Gross-Neveu-type is: 

2 det A A12 

'^^"l6(l + ,^)A22ln(Ai)' '^^^ A22' 



(59) 



(60) 



The solution of the mass-gap equation of the Abnormal-type is: 

4 . ,0 A^/'(3y3Ai2 - A22)'/=^ 



9(1 + ^7^) -^4. 35/3(1 + ^4) 1^2/3 (Ai) 



(61) 
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For the case 2, of a non-zero p , the mass-gap equations reads: 

A2 



All 



Ai2 

A22 



160^1n-_ + 16(7V?ln 

'^1 



+ (1 + ^?^ 



159 



5^3 



'102 + 7 (02 + 



5v3 ,2 3 , , v3 / ,2 

■^0? + 70102 + ^ (02 + 



9 



1 + 9 ) T0I + ^0102 + -Ar2 + p 



(62) 
(63) 
(64) 



The mass-gap equations (|6^-(|64D can be rewritten in an equivalent form: 



A2 



3An - 3A22 + 9 V3A12 = 480? In —2 - 3 + 48^V? In 



A^ 



A22-3V3A12 = (1 + (?^)(1202- 4^30102), 

^0? + ^0102 + ^(02+P' 



^22 



[1 + 9' 



(65) 
(66) 
(67) 



From (|65|)-(|67D it is clear that for fixed A^/ while A grows larger, the solution exists if 
Afc/ parameters are chosen close to a particular plane in the parametric space. This plane 
is defined by the equation: 

A22 = 3v^Ai2. (68) 

When Afc/ satisfy the equation ( |68|) exactly, the solution is found to be as follows (in the 
large-log approximation): 



An 



16(1 + g^) In 



/,2 



1 + 



In 



A2 



4 A 



22 



9(1 + ^7^ 



V301, 



1 + 



In^ 



(69) 

(70) 
(71) 



4.2. Mass spectrum in Model I 

The mass spectrum of related bosonic states is determined by the Eqs. (26)- (31) and 
taking into account the conditions necessary for a minimum of the potential (52,53). The 
solutions at — = < one can obtain from the Eqs: 

det(Ap2 + 5) =0, (72) 

The "kinetic" matrix A as being proportional to is derived in the soft-momentum 
expansion in powers of and in a large-logA approximation. Because the expressions 
for A and B are cumbersome we give their explicit form in the Appendix A and B 



16 



correspondingly. After substituting expressions for the matrix A, B into (^) one can get 
the mass spectrum for the neutral Higgs bosons in Model I. For the case, p = 0, the mass 
spectrum is resembling ones in Two-Channel model, in particular, the Gross-Neveu-type 
solution brings the spectrum for scalars: 



det A 



~ 77IY = (74) 

(l + ^?2)A22ln(f) 



and for pseudoscalars: 



ml = 0. (76) 
The Abnormal solution induces the following mass spectrum for scalars: 

« (77) 



and for pseudoscalars: 



~ , / = 6m^ (78) 

32/3(1 + (^2) 1^2/3 (^A^^ ^ ^ 



^ 32/3(3V3A,,-A,,)V3 

27(1 + , 2) aV3 1,1/3 (Ai^' ^'^^ 



m2 = (80) 

Remark that dynamical mass irid is in fact the mass of t quark in the Model-I, because 
the v.e.v. of which is parametrized as 0i = gives the value of mass of t-quark. 

The mass spectrum in the P-parity Breaking Phase, for a non-zero p is: 

ml = 0, (81) 



ml ^ ^ ^ 160? = 4m2, (82) 

' 3(l + ,2)inA^ v^i V ; 

fj, 

2 ^ SA22 , X 



3(1 + ^2) 

2 3An + 7A22 

" 27(1 + ,2) In Ai' 

Thus, we have constructed the Model I where: 



a) Two composite Higgs doublets are created dynamically as a consequence of DCSB 
in two channels. 
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b) In 2HQ Model I Higgs bosons are rather radial, ground and excited states in the 
scalar-pseudoscalar channels. 



c) The appropriate fine tuning leads also to spontaneous breaking of P-parity and, 
therefore, of CP-parity in the Higgs sector. 



5. Top-Bottom Condensation for 2HQM Model- II 



5.1. Effective potential in Model II 

The Lagrangian density of the Model II to describe the dynamic of two composite 
Higgs bosons which consist of bound states (condensates it,bb) and satisfy the FCNC 



suppression |]T6| can be written as: 



2 2 



^5) 



k,l=l 



where Pl{R) = 1/2(1 ±75) - the left and right projectors, and M is the two-by-two fiavour 
matrix: 



m=l 



v 



^m2 ft,m I "^^2 j 0ml ft,m I "^^2 



-'ml fb,m 



A2 



0m2 fb,r. 



Pi] 

J 



(86) 



where we set for the Yukawa coupling constants gt,k,gb,k = 1 (two Yukawa constants 
are due to the renormalization of Higgs fields and other ones we choose equal one). In 
this Model II $i,$2 give masses to up— , down— type quarks. The structure of quark 
interaction is specified in four formfactors: 



/t,i 

ft,2 
fb,l 
fb,2 



1 - 



^7) 



When the chiral symmetry is broken, the v.e.v. of neutral Higgs fields are non-zero 
and the true Yukawa vertices should be obtained by subtracting from M its v.e.v. 



M = M-(M), 



where (. . .) means a v.e.v.: 



(M) 



TTlt 

e'^°mb 



^9) 
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The elements of the matrix (|89|) are the quark mass functions: 

mt{T) = 0i(l + Q,ir) + 02e*'^°Q,2r, 

m5(r) = 0iCfe,ir + 02e*'5o(l + Cfe,2r), ^ = -J2^ (9°) 

defined to be real and 0i =< 0i2 >,02 =< ^22 >• The non-zero phase at rrib, which is 
displayed explicitly in (^9|), may appear if the v.e.v. of ^22 acquires irremovable phase 
factor when the chiral symmetry is broken. 

As the vacuum charge stability is assumed, (M) is diagonal, so (M) and (M)''" commute 
and can be placed in any order in products of themselves. 

The effective potential of composite two-Higgs model II in which the interaction of 
quarks and Higgs bosons is described by formfactors (|^) reads: 



1 / 1\ 

Veff = -An0?-2Ai2(/'i(^2Cos5o- A2202 + 2*^1 + 2 1 + 

1 / A^ l\ 1 
+ -0^ Mn^ + -j + -(Jiiii0f + 4Jiii20?02Cos5o + 

+ 2(Jii22 + J122I + -^1212 COs25o)0i02 + 4^12220102 COS 5o + 

+ J22220^)+O(^), (91) 

where v.e.v. of fields are: < 0i2 >= 0i, < 022 >= 02- (For the concrete choice of 
formfactors in Model II the form of the effective potential is displayed in the Appendix 
C) Its minimum is described by solutions of the mass-gap equations (21), (22) which for 
the Model II are: 

2 Aii0i + 2 Ai202 cos 5o = 

= 203 In 1^^^ + 2Jiiii0? + 6Jiii20?02 COS(5o + 2Ji222 02 COS(5o + 

+ 2 ( Jll22 + -^1221 + -^1212 COS 2So) 0102, (92) 

2Ai20i02sin5o = 

= 2 Jiii20i02 sin 6o + 2 Ji2i20i02 sin 26o + Ji2220i02 sin 60, (93) 

2 A2202 + 2 Ai201 COS 5o = 

= 202 In ( 02 j + 2-^2222 02 + 6 Jl22201 02 COS + 2 Jiii20i COS + 

+ 2 ( Jll22 + -^1221 + -^1212 COS 25o) 0i02, (94) 

where Jkimn {k, l,m,n = 1, 2) are the integrals: 

1 



+ fb,k (r) fb,l {T)fb,m{r)fb,n{T) + 
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+ft,k (r ) fb,i [r) fb,m (r) /t,„ (r) - 

- /t,fc (t) /t Am (t) /fe,n (t) + 

+fb,k (r) (r ) /i,^ (r) (r) - 

- /b.fe (t) ft,i (r) /t,m (r ) (r ) - 
-Afc(0)/M(0)A™(0)/i,„(0)- 

-/t,fc(0)/M(0)/fe,^(0)/t,„(0) + 

+Afc(0)/M(0)/fe,m(0)/M(0)- 
-/6,fe(0)/M(0)Am(0)/t,„(0) + 

+Afc(0)/u(0)A„(0)/b,„(0)) -. (95) 

r 



It is more convenient to solve the equations for the variables Aim rather 

than 01, 02, ^0- The variables 0i, 02, So will be treated as input parameters while Aim 
as the unknowns. The reason for this is that we do not know A^^ from any global theory, 
we rather fit them so that 0i, 02, Sq conform to experiment. The equations (|92D-(p4D 
are linear for Aim and can easily be solved; one just expresses Aim. through 0i, 02, 6o 
and substitutes them in every place they appear. As usual two cases must be considered 
separately. 

1) For 5o = 0: 

All = 0? In (^^^ + Jllll0i + 3 Jiii20102 + 

0"^ 02 

+ (^1122 + <^1221 + <^1212) 02 + ^1222" ^12 — , (96) 

01 01 

A22 = 02 In + ^2222 02 + 3 Jl22201 02 + 

+ (J1122 + J1221 + J1212) 0? + ^1112^ - Ai2^, (97) 

02 02 

where A12, 0i, 02 are treated as input parameters; 

2) and for Sq ^ 0: 

All = 0? In + Jiiii0? + 2 Jiii20i02 cos 60 + 

+ (-^^1122 + -^^1221 — -^^1212) 02, (98) 
A22 = 02 In + J2222 02 + 2 Jl22201 02 COS Sq + 

+ (^1122 + ^1221-^1212)0?, (99) 

A12 = ^111201 + -^1222 02 + 2 Jl2120102 COS (5o. (100) 

The mass spectrum of related bosonic states is determined by the matrices A and B of 
the second variations of the effective potential ( ^) (see Appendix A, B). 
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5.2. Mass spectrum in Model II 

After substituting explicit forms for the A,B into (|72D, one can obtain the mass- 
spectrum for the composite neutral Higgses in Model II. 
1) For 5o = 0: 





^ 2mt, 




= 0, 











Ai2 - J1112) 



(101) 
where 

r=!I!lln^^l. (102) 
If the ratio were (say, for the fourth generation): 

777 

^ = 0(1) (A ^00), (103) 
rut 

than one would get: 



2 _ 2(mg + ml) 
mfTrif^ in 



2) For 60 ^ 0: 



X A12 - J1112 - 2mjm^Ji2i2 - 'm^Ji222 ■ (104) 



(105) 
(106) 



mi 


= 0, 


m2 


K, 2m 


ma 


^ m^ 



^ (107) 
\ In 4 

m4 2rm J sin 5o 1^71212 . (108) 

For the case of 5^ 7^ 0, the model predicts a low value for 171^. 

If one considered (that may take place for the fourth generation): 

~ m^ ~ 1 (A ^ 00), (109) 

the mass-spectrum would turn out to be as follows: 

(110) 
(111) 
(112) 



mi 


= 


m2 


2m^ 


ms 




m4 


2 sin (5 



Juu{m^ + ml)_ (113) 



In 

mt 
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We notice that when Sq we have neither scalars no pseudoscalars any longer because 
the particles which are eigenstates of the energy operator, are mixed of both P-even parity 
and P-odd parity fields, hence the former classification by parity does not hold for this 
particular case. 

5.3. Mass spectrum in Model II for the special choice of Yukawa 
couplings with FCNC suppression 

The suppression of the flavor changing neutral currents in a composite two-Higgs model 
can be implemented when Yukawa couplings are chosen as follows: 

gt,i = gb,2] gt,2 = gb,i = \/5. (114) 

It can be easily seen from the potential in Appendix C that with the above choice of 
gt,k, gb,k the last two terms in (191) come to zero. This prevents the appcarcnce of the 
Yukawa coupling inducing FCNC from radiative corrections. In this case the potential of 
Model type II reads: 



eff 



Aii($l$i) - Ai2($l$2 

A2 



87r2 ^ 



A2l($^$l) - A22($^$2) + 



. v+(0) 2 2 



+ 



33 3 

+ -($l$l)($^$2) - -($I$2)($^$l) + 



+ - 



15 



;$I$2)' + (*I$2)' 



:ii5) 



where for i^±(0) we have 

t/± = 2($l$i) + 2($^$2) ±2[{^\^if 4 
±2($l$i)($^$2) - 4($l$2)($^*i) 



1/2 



(116) 



To be short we skip the calculations and display only the mass spectrum. When the 
mass ratio of t- and 6-quarks rrit/mi) holds unchanged while A goes to infinity, in the 
large-log approach we get the following estimations for Higgs' masses: 

1. for the case of concerved P-parity {5q — 0): 



= 2mb, 
^ 2mt, 
= 0, 



+ ml) 



2m*mft In ^ 



15 

T 



Ai2 - —mtm^ 



(117) 
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2. and for the P-parity breaking phase: 



mi 


= 


m2 


~ 2m^ 


ms 


^ 2mb 


m4 


sin (5 



\ 



15 (mf + m^) 



8 In 



A2 



(118) 
(119) 
(120) 

(121) 



Thus in the Model II the Quasilocal Yukawa interaction with Higgs doublets reduces 
at low energies to a conventional local one where each Higgs doublet couples to a definite 
charge current and its v.e.v. brings the mass either to up- or to down- components of 
fermion doublets. Based on the FCNC suppression, the Model II leads to the relation 
rrit » rrih and so to an enhanced coupling of the light scalar (pseudoscalar) boson to 
the down-type quarks while suppressing the coupling to the up-tjpe quarks. The Model 
II has a broad spectrum of excited bound states which is to parametrize the data, in 
particular, obtained from the Next Linear Collider. 



6. Summary 

In our paper we have proposed a set of Quasilocal NJL-type quark models (QNJLM) which 
lead to a larger spectrum of ground and excited states in the polycritical regime. From the 
viewpoint of the SM, these models are considered as more natural than common extensions 
of the SM, since they do not enlarge the number of elementary particles in fermionic sector 
and preserve the symmetries of the SM. For the toy Two-Channel Quasilocal quark model, 
near tricritical point we have found three major phases: a symmetrical one and two phases 
with DCSB, different in correlation lengths in scalar channels. On a particular plane in 
the space of coupling constants we discovered the special P-parity breaking phase. It 
means that in such a phase there exist heavy scalar states which can decay into two 
or three light pseudoscalars. This phenomenon of dynamical P-parity breaking can be 
used in the extensions of the SM where several Higgs bosons are composite ones. In the 
framework of the QNJLM we have presented two Models which provide at low energies 
two composite Higgs doublets, as minimal extensions of the Top-Mode Standard Model 



ri|],[jT4|. In the 2HQ Model I Higgs bosons are rather radial, ground and excited states 
in the scalar -pseudoscalar channels. In the 2HQ Model II, which consistent with the 
requirement of natural flavour conservation strong forces lead to the formation of 



top and bottom bound states (and corresponding condensates) and generate masses of 
t,6-quarks. In Model II we have concentrated on the scenario where each of the neutral 
components of the two doublets 0i^2 (with v.e.v. fi 2) couple at low energies respectively 
to the /a = ±1 fermion flelds The FCNC suppression leads to the relation nit » nib 
and to an enhanced coupling of the light scalar (pseudoscalar) boson to the down-type 
quarks and the charged leptons while suppressing the coupling to the up-type quarks. The 
existence of light neutral Higgs (pseudo) scalar bosons in the framework of 2HDM is not 
excluded by existing data (< 40 GeV). The chance that it can be seen at the Next Linear 
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Collider in the 77 processes has been pointed out in |g6|,|]2^. As a result of complexity of 
two v.e.v.'s for two composite Higgs doublets the dynamical CP- violation may appear in 
the Higgs sector. At high energies these channels are strongly coupled and one could say 
that two-composite Higgs doublets partially represent the mixture with excited states. If 
such excited states exist then they will modify the Higgs mass predictions. In addition, 
we remark that low values for the Higgs masses of the additional states could actually 
change the window for since these states could give a significant contribution to the 



p-parameter p8[. From our consideration we have seen that the appearance of dynamical 
CP- violation in the Higgs sector imposes strong bounds on Higgs masses, in particular, 
one light scalar Higgs boson is unavoidable. The experimental implications of such effects 
are expected to be rather small in the fermion sector of the SM These effects are 

observable in decays of heavy Higgs particles (namely, pseudoscalar Higgses may decay 
into scalar ones, scalar Higgs may decay into pseudoscalar ones) and in decays of Higgses 
particles into two vector bosons where CP-even and CP-odd amplitudes appear. At 
high energies the appearance of the appreciable CP-violation could be important both 
as a source of electron and neutron electric dipole moments and as a mechanism for 
EW scale baryogenesis []30[| , []3T| . Besides one expect also that modifications of the SM 
Lagrangian (the Higgs and Top interactions) by higher dimensional vertices may enhance 
the Higgs production at hadron colliders |]32| . 

The theory of two composite Higgs bosons which we have discussed in our paper should 
be regarded as a viable alternative to other approaches to the BSM and perhaps the major 
progress in the alternative approaches will come when the first direct experimental results 
associated with the origin of EWSB begin to appear. 

The purpose of this paper has been to elaborate the very design of quasilocal NJL- 
quark models with two-composite Higgs bosons. A more comprehensive analysis of low- 
energy particle characteristics in these models is postponed to the next paper in this series 
of. The numerical computation of bounds on mass spectra, Yukawa coupling constants 
and decay widths with taking into account the renormalization-group corrections will be 
presented elsewhere. 
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Appendix A: Kinetic matrix A for composite two-Higgs 
bosons 

In this appendix we calculate the kinetic term for composite Two-Higgs Quasilocal 
Quark Models which is obtained by calculation of the one-loop diagram: 




Figure 1: One- loop diagram for calculation of kinetic term. 

Here p is an incoming momentum, and A; is a momentum running around the loop. 
The loop diagram (Fig.l) gives the following expression: 



|fc|<A 



(2^ 



Tr 



(iMPL + iM^Pn) D"^ [ (k + 



:i22) 



Here the coefficient 1/2 is due to symmetry of the diagram (Fig.l). The full expression 
within the square brackets in ( |122| ) is an element of a direct product of three spaces: color, 
flavor and spinor and the trace operation is related to all them. We define the vertex: 



z{MPl + M^Pr) 



and the fermion propagator: 



(g-z(M)Pi-2(M)tp^)A(g2), 



where A is a matrix function: 



A(g2) = (g2 + (M)(M) 



2 2 / ^ 



1 



(123) 
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First, we calculate the trace of the sum of all the products of 7- matrices displayed in the 



expression (|122|) . After that we come to the following expression: 



2 



|fc|<A 

-2tr 
+2tr 
-2tr 



(2^ 



2tr 



MA( (A: + | 



M'Al(k-l 



+ 



(124) 



M(M)tA ( ( A: + I j J M(M)tA [(k - 



MA 

2^ 



P 



P 



In this formula and further on the trace is applied only for the flavour two-by-two matrices. 

The kinetic term is derived as being proportional to in the soft-momentum expan- 
sion of (|124|) in powers of p'^. We obtain this term by means of calculating the second 
derivative of ( |124| ) at zero external momentum p. First, let us rewrite the expression ( |124| ) 
in a form: 



_Nc 
327r4 



fix,y) - phi 



where 



|fc|<A 



X 



mA{k'^)yi A{k'^ 



k + 



P 



k- 



P 



(125) 



(126) 



with the function / defined as: 

f{x,y) = 4tr 

- 2tr 

- 2tr 



MA(x)M A(?/) 



M(M)^A(a;)M(M)^A(y)" 
m\m)A(x)m"^(M)A(?/) 



(127) 



Let us expand the expression ( |127|) in series of p and extract the term proportional to p"^ : 



1 92 



32n'^ J \2 dpodpo 



fix,y) -tr 



p=0 



MA(A;2)M A(A;2 



(128) 



where the angular brackets stand for angular average in 4-dimensional Euclidean space 
and Pq is a component of 4-momentum p. The second derivative of the function / reads: 



, dpodpo 
1 
2 

+fx{xo,yo) + fyixo,yo) 



p=0 



{fxx{xo, yo) - 2/^j,(xo, yo) + fyy{xo, yo)) y 



(129) 
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The subscripts in fxx, fxy, fyy, fx, fy stand for partial derivatives by variables x and y. 
The derivatives are calculated at xq = yo = k"^ (see (|126|) ). The function A, defined in 
( |123| ), is a flavour matrix: 



/ 



A(x) 



1 







X + m. 



A2 



X + 



Using m 



f{x,y) 



I), O and (|T30I), one gets 

/ ft,l (az) ft,m (a^) 



E 



Lm=l K - 



_ {x + mt(P/A2))(?/ + m;,(P/A2)) 



+ 



+ 



fb,l (f^) /fc.r 



X + mfe(A;VA2))(y + mt(A;VA2))y 



A2 

2 



A2 i 



-i5o 



A2 



+ 



+ h,i 



p_ 

A2 



t,m 



A2 



+ 



+ 4k 



+ 



^ {x + mtikyA^)){y + mb{kyA^)) 

I ) 

X + m6(P/A2))(7/ + m* (P/A2)) J 

fb,l (a^) fb,m (a^) 



+ 



(x + mfe(A;2/A2))(y + mb(fc2/A2)) 

ft,l (aS') ft,m (a^') 



+ 



+ 



{x + mt(fc2/A2))(y + mt(P/A2)) 



- 2 



2i<5o 



/m (^) (g) (1^) e-^- 
(x + mfe(P/A2))(y + m6(P/A2)) 

ft, {^) km (^) 
X + mt(A;2/A2))(y + mt{k'^ / K^)) 

fb,i (^) (f^) m2 (Ai) e2^^o 
(x + mfe(A;2/A2))(y + mfe(P/A2)) 

ft,l (fj) ft.m (a2) "^t (^2 



+ 



+ 



(130) 



(131) 



(x + mt(P/A2))(?/ + mt(P/A2))_ 
After applying the derivative procedure displayed in ( |129| ) to the function /, one gets the 
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kinetic term in the following form: 



N ^ / 



16^' l,m=l V 

+l\^{d,ct>n){d,ct>m2) + /S(5.0r2)(9A)) ■ (132) 

contributes to the kinetic term for charged components of higgs doublets, //^\ and 
do the same for the neutral components. The expressions for them are cumbersome 

and we have divided the total expression in three parts. For the charged components one 

has: 



^Im — 



I- ^ 



^2 ' 



mi{kyA'')e 

^ + m?(A;VA2))(A;2 + ml{k''/K^)Y ^ 

k^ 

^ (A;2 + m^t{kyA^)y{k^ + ml{k''/A?)y ^ 

+ \ U 

(P + m2(A;VA2))(A;2 + m2(A;2/A2))y 

X [f't^l (f^) ft,m [j^ + h,l (^1^ j h,m (x^J J + 

^ V(A;2 + m?(fc2/A2))3(fc2 + ml{k^/A^)) ^ 

mf(fcVA=^)m,(A;VA^) 
^ (F + m?(A;2/A2))(A;2 + ml{k''/A^)f ^ 

m,{kyA^)mt,{kyA'')e \ 
^ + m?(A;2/A2))2(A;2 + m2(A;2/A2))2 J ^ 



fc^d/c^ (133) 



for the neutral ones: 



/■(2) ^ 1 



A2 



Jt,l J yj^2 J y (^2 + ^2(fc2/A2))4 ^ 



2fc2 1 \ 

^ (A;2 + m|(A;2/A2))3 + (A;2 + m^t{kyA^)y ) ^ 
+ A' (x^) A- (x^) (~(A;2 + m2(A;2/A2))4 + 
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+ 



+ 



and 



7-(3) 



1 "^^r / 
2/ ^*''( 



A2 j J''"^ \^A2 y \^2(A;2 + m2(A;2/A2))4 + mt2(fc2/A2))3^ 



A:2m2(A;2/A2) 



m2(A;2/A2) 



(134) 



+ 



2i(5o 



2(A;2 + m2(A;2/A2))4 



m2(A;2/A2) 



(A;2 + rng(A;2/A2))3 



(135) 



7-(4) 



1 ^'r / 

2/ ^*''( 
n L N 



/ fc2m2(fc2/A2) m2(fcVA2) 

A2 j |^A2 J I^2(fc2 + rnf (fc2/A2))4 + mlik-^/K')) 



+ 



m2(A;2/A2) - 
(F + m2(A;2/A2))3^ 



2(fc2 + m2(A;2/A2))4 



(136) 



The next task to do is to calculate the integrals for the large value of A, ignoring all 
contributions which disappear in the A — > cxo limit. Thus one obtains: 



Ii^ = (/t,K0)/t,m(0) + /M(0)/6,m(0))x 



X 



A2 Smlm^ — mf, , m'} mi + Grn^ml + mf 



In h 

rrit {rrit — m2)3 



In^ 

mi 



A{m'i — mlY 
+ (/M(0)/b,„(0)e-^« + /^,K0)A„^(0)e*'5«) x 



+ 



X 



2mlml ^ ml m^mf^i^m"^ + mf) 
ml A{ml — mlY 



[ml — mIY 



+ 



+ j {ft,l{^)ft,m{r) + fb,li^)fb,mir) - 



dr I In 



A2 



(137) 



= /M(0)/t,m(0) 



In 



A^_ 13 
m? 12 



+ AKO)Am(0) 



In 



A2 



13 
12 



+ 



+ 



j Ut,i{^)ft,m{r) + fb,i{r)fb,m{r) 



(138) 
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ilf = -liMo)fUo) + hMhA^y-^''') + o I 1 , (139) 




4? = -^(/t/o)/,„(o) + hMhA^y'-'') + oy^\. (140) 

Herein and further on mj and mj, stand for quark masses. 

As we know from experiment, the mass of the top quark is much greater than that of 
the bottom quark. Regarding ^ m^, for the choice of formfactors (|87D one gets: 

^ \n^-l + 2ct,i + hcl + cl,), (141) 



Iv2 ~ Ct,2 + C6,i + -(Ct,iQ,2 + Cfe,iCfe,2), (142) 

/« ~ inA!_^ + 2c,,2 + ^(42 + ^2), (143) 
mr 4 2 ' ' 



111 ^ l^^-i + 2c*,i + ^«i + CM), (144) 

(2) 1 

112 ~ Ci,2 + Cfe,! + -(Ct,iQ,2 + Cfe,iCb,2), (145) 

12 2 



1:^2' ^ In— - — + 2c6,2 + -(c'2 + <2), (146) 



lif ^ -I, (147) 

lif ^ 0, (148) 

/f «^ -^e-2^^°. (149) 

4? - (150) 

4? ^ 0, (151) 

4? - -^e^'^o (152) 

Next, we shall change variables and rewrite the total expression for kinetic term. For the 
neutral components we choose non-linear parametrization: 

0ie^° = 012, 026*^"+'^ = 022. (153) 



Let us substitute (|153| ) into ( |132| ) and expand it with derivatives, leaving only quadratic 



terms and omitting the rest of the expression. As the v.e.v. of 0i and 02 provide quarks 
with dynamical masses, we replace them with nit and rrib respectively. The variable a is 
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regarded as a Goldstone boson, which is absent in the effective potential; it appears only 
in higher-derivative terms. The other phase, S , is associated with the relative phase. The 
variables 0i and 4>2 parameterize radial excitations. For the fields a and S we use different 
notations 

5 = 03, a = 04, (154) 

so that one can rewrite the kinetic term for the neutral components of Higgs doublets 
uniformly: 



J2 Am{d^,(j)l){dij,(j)m), 



where A is four-by-four matrix: 



l,m=l 



(155) 



.(2) 2 

"3 


lif cos So 


-lifrrib sin So 


7-(2) • X 

-/iVmfe smdo 


lif COS So 


7-{2) ^ 
-'22 ~ o 





42^ sin So 


7-(2) • X 

- mfosm do 







(4? + 

(2) 

-|- rrihiTi^ cos^o 


r(2) • X 


lil^rrit sin So 


(2) 

lu m^m^ cos 5o+ 
+ (4? + 


(4? + |)m?+ 

-|- 27^2^ m^m^ cos (5o+ 

+(4? + i)m^ 



(156) 



Appendix B: Momentum independent matrix B 

Let us define the matrix of second variations of the effective potential for the Model I in 
the following way: 

For the case, when p = 0, the matrix for scalars and BIJ^ matrix for pseudoscalars 
are represented: 



^11 



(1+/) 



A2 



15\/3 



Vs. 



+ 



TDcra 
^12 



^22 



^11 



01 



-2Ai2 + (1 + g^) 
(1 + 



902 + 

'5^3 , 



90102 - 
0102 + 



2 2 02 



15V3 ,2 ^ 3V3 2 

^01 + ^<^2 

5^0? 



+ 2A 



12 



o^2 



^<^102 ~ ~ 302 



02' 
01 



(158) 
(159) 
(160) 
(161) 
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DTTTT 
-D22 



^kl 



= -2A 



12 



"5^3 ,2 



^1(^2 



A;,/ = (l,2). 



02 



(162) 

(163) 
(164) 



For the case, p non-zero, the corresponding matrix of second variations of the effective 
potental is a 4 x 4 matrix. One can arrange the neutral Higgses in a vector-column with 
4 components: 

0-2 

where Ufc - scalar fields, tt^ - pseudoscalar. 

Let us display the non- vanishing components of the matrix of second variations, cal- 
culated in the minimum of the potential, in the following form: 



ID a a 
^11 



Tjaa 
-D12 

Tjaa 
-D22 

-D12 
-D22 

JDCriT 

-D12 

rjcrir 
^22 



640^ In 



- 22302 - 10^301^ 



-640?^Mn 

= (1 + /) [-5V^0? + 60102 + ^0: 

= (1 + /) [30? + 2v^0i02 + 90: " 

= 3(1 + <7V, 

= v^(l + /)p^ 

= 9(l + /)p^ 

= p(H-(7^)(-5x/30i + 302), 

= p(l+/)(30i + V^02), 

= p(l + (?^)(30i + V302), 

= p(l + ^^)(v^0i + 902). 



(165) 
(166) 

(167) 

(168) 
(169) 
(170) 
(171) 

(172) 

(173) 

(174) 



The matrix B for the effective potential in the Model II is: 
Bu = -2An + 60? In ( ^ ) - 40? + 6 Jiiii0? + 



^1, 

/,2 



-M2 Jiii20l02 COS 60 + 2 ( J1122 + -^1221 + -^1212 COS 25o) 02, (175) 

i?i2 = -2 Ai2 cos 6o + 6 Jiii20? cos 5o + 6 Ji222 02 cos 6o + 

+4 ( Jii22 + Jl221 + -^1212 cos 26o) 0i02, (176) 

Bis = 02 sin 6o (2A12 - 6 Jiii20? - 8 Ji2i20i02 cos 60 - 2 J122202) , (177) 

Bu = 0, (178) 
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B21 = Si2, (179) 

B22 = -2A22 + Q(pl In l^li j - + 12 Ji2220102 COS Sq + 

+6 J2222(/'2 + 2 (^1122 + ^1221 + ^1212 COS 26o) 0i, (180) 

B23 = 01 sin 5o (2 A12 - 6 Ji22202 - 8-^12120102 COS Sq - 2 Jiii20?) , (181) 

S24 = 0, (182) 

Bsi = Si3 (183) 

^32 = S23, (184) 

-B33 = 20102 (A12 COS So - Jlll20i COS Sq - 2 Ji2l20102 COS 2So - 

-^1222 02 COS (^o), (185) 

B41 = 0, (186) 

542 = 0, (187) 

^43 = 0, (188) 

S44 - 0. (189) 



Appendix C: The effective potential of Model II for 
the special choice of formfactors 

For purposes of further calculations of realistic mass spectra, Yukawa coupling constants 
and decay widths with taking into account the renormalization-group corrections we 
present in this appendix the effective potential of Model II for the following set of form- 
factors: 





/<,i = 2-3r; f,^, = -V^r; 






fb,i = -\/3t; fb,2 = 2 - 3t; 


(190) 


the constants Jkimn 


are evaluated to definite numbers. Seven of then are defined as follows: 


■-^1111 — 


^(341-53^1), 


(191) 


■^1112 = 


^ (^6,1^6,2 - 5gl^9t,2) , 


(192) 


■-^1122 — 


I {glA - 53^6,2^4,1 - '^9b,i9b,2gt,igt,2 + ^gli9l2 + qIiqI'^ , 


(193) 


•^1212 = 


3 / 2 2 2 2 \ 

^ \9b,i9b,2 + 9t,i9t,2 ) ) 


(194) 


-^1221 = 


1 (^6,1^6,2 + 53^b,2^t,i + '^9b,i9b,29t,i9t,2 - ^9li9t,2 + 9t,i9t,2) , 


(195) 


-^1222 — 


^ {9t,i9t,2 - 556,156,2) , 


(196) 


<^2222 = 


1 (35^2 - 53^j2) • 


(197) 
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The rest of them is found from their symmetry property: 

Jklmn Jmnkl Jlknm- 

With ([T9TD-([T9^) the potential for the Model II reads: 



-An($I$l) - Ai2($I<I>2) - A2l(<I>5$l) - A22($5$2) + 



In 



^fln4, + lU 
2 I u_{0) 2 



+1 {Hi - 53gt) + I (3<2 - 53^7,^,2) (^1^2)' + 



+ 4 (^?M^?6,2 - 53(7',2^7m - 2(?b,i(?fc,2(7t,i(7t,2 + 39I19I2 + olA) ($I$i)($|$2) + 
+ T fe,i^6,2 + 53^fe%^2i + 2gbAgb,29t,i9t,2 - 39li9t,2 + ^t,i^'2) ($I$2)($5$i) + 



4 
3 

'4 
/3 

+ — {9!,i9b,2 - 5gl^9t,2) [($i$2) + 



+ 



/3 



(19^ 



where we adopt the definition for i'±{0): 

= 2^7ji($l$i) + 2^7,%($^$2) ± 2[gli<!>\<^^f + ^7j2('^'2<^'2)' + 
+2glgl,i^{^,)i^l^2) - AglglM^2)i<^l<^i'''^' 



(199) 
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